import statsmodels.formula.api as smf
import bambi as bmb
# Fit a frequentist logistic regression
fr_fit = smf.glm(
'y ~ x1 + x2 + x3',
data = conjoint_data,
family = sm.families.Binomial()
).fit()
# Fit a Bayesian logistic regression
ba_fit = bmb.Model(
'y ~ x1 + x2 + x3',
data = conjoint_data,
family = 'bernoulli'
).fit()

